Pressure-induced phase transition and bi-polaronic sliding in a hole-doped CU2O3 
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We study a hole-doped two-leg ladder system including metal ions, oxygen, and electron-lattice 
interaction, as a model for Sri4_ I Ca a; Cu24 04i_5. Single- and bi-polaronic states at |-hole doping 
are modeled as functions of pressure by applying an unrestricted Hartree-Fock approximation to a 
multiband Peierls-Hubbard Hamiltonian. We find evidence for a pressure-induced phase transition 
between single-polaron and bi-polaron lattices. The electronic and phononic excitations in those 
states, including distinctive local lattice vibrational modes, are calculated by means of a direct-space 
Random Phase approximation. Finally, as a function of pressure, we identify a transition between 
site- and bond-centered bi-polarons, accompanied by a soft mode and a low-energy charge-sliding 
mode. We suggest comparisons with available experimented data. 
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I. INTRODUCTION 



A material of considerable recent interest is the com- 
pound Sri4_ :E Ca a ;Cu2404i (SCCO), which shows super- 
conductivity under pressure 0,0 The SCCO structure 
includes quasi-one dimensional (Q1D) two-leg ladders 
CU2O3 and one-dimensional (ID) Cu02 chains, while 
other high-T c cuprate superconductors found so far con- 
tain two-dimensional (2D) CuC>2 planes. The CU2O3 
ladders and Cu02 chains in this system are intrinsi- 
cally hole-doped even at x = with a total of 6 holes 
per formula unit. (The total hole concentration can be 
decreased by La and Y substitution for Sr, with e.g. 
LagCu804i containing no holes.) Of these, only approx- 
imately one hole goes into the ladder component of the 
formula unit, (0^03)7, which results in an effective dop- 
ing of about 7% per Cu site in the ladder (3, 0. Ca 
substitution, x, does not change the total number of car- 
riers, but transfers the holes from the chains to the lad- 
ders 0, 0, @ • The conductivity increases with increasing 
10. 

The doped holes can create a polaronic or charge- 
density-wave (CDW) state, and a charge sliding mode 
could be expected as a collective excitation. The exis- 
tence of such states and modes is supported by some 
experiments in the material with x = 0. Resonant X-ray 
scattering has revealed a five-site periodic hole structure 
in the ladder |8|, |9( • Microwave measurements show a rel- 
atively small c-axis conductivity with a narrow peak in a 
very low-energy region (~ 0.2 meV) OS El- 

This low- 
energy resonance is observed up to a temperature (~ 100 
K) too high to be attributed to single particle excitations, 
which would be completely broadened and no longer ob- 
served above 10K (~lmeV) due to thermal fluctuations. 
Similarity has been noted between the nonlinear behav- 
ior of the conductivity in SCCO at x — and that of the 
sliding mode in materials s upp orting CDW states. Blum- 
berg et al. have reported [l^] that the low-frequency di- 
electric constant (eo ~ 4 x 10 6 ) obtained by transport 
measurements is consistent with estimates from the pin- 
ning energy (~ 0.2 meV) suggested by microwave mea- 



surements. 

There are several possible polaronic configurations. 
One is to randomly distribute the polarons. Another is 
that some of them bind to compose various multi-polaron 
configurations: In a 2D case, the doped holes can arrange 
into stripes, which is one kind of multi-polaron state. The 
groundstate configuration depends on the doping level, 
pressure, etc. Modeling the pressure-dependence of the 
groundstate configuration is one aim of this study. 

In a previous work 0, 0] , we studied the electronic 
and phononic excitations in a 2D Cu02 plane with in- 
homogeneous charge-lattice-spin structures (stripes and 
other polaron patterns). We identified local "edge" or 
"interface" modes in the phononic and electronic (spin 
and charge) excitations induced by the inhomogeneity. 
In the CU2O3 two-leg ladder in SCCO, we can similarly 
anticipate the existence of inhomogeneous structure and 
associated local excitation modes in spin, charge, and 
lattice degrees of freedom. 

Here, we investigate the ground states and the elec- 
tronic and phononic excitations in the two-leg ladder 
system by applying an unrestricted Hartree-Fock and a 
direct-space Random Phase approximations (RPA) to a 
multi-band Peierls-Hubbard Hamiltonian. We consider 
single-polaron (SP) and bi-polaron (BP) states: The lat- 
ter comprise bound polarons extending over rungs. Com- 
pared to the 2D cases, the SP state that includes iso- 
lated polarons, is found to possess similar phonon excita- 
tions as the diagonal stripe state or the periodic polaron 
state 10, 0| . On the other hand, the BP state shows the 
same type of local phonon modes as the vertical stripe 
state in the 2D system. This is reasonable, since the ver- 
tical stripe state is a form of multi-polaron state, which 
includes several polarons coupled by shared O ions. To 
model the effect of pressure on the groundstate configu- 
ration, we compare the energies of the SP state and the 
BP state, while varying the Cu-0 hopping integral. We 
find that as a function of increasing "pressure" (modeled 
by increasing hopping strength) a transition from the SP 
state to the BP state is induced, together with interesting 
intermediate states. Most strikingly, we identify a tran- 
sition between site- and bond-centered BPs accompanied 
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by phonon softening indicative of the onset of sliding or 
other instabilities. 



the phonon eigen-modes, we calculate the corresponding 
neutron scattering cross section: 



II. FORMULATION 



dt e~ luJt 



ikrf(0) e ikr, e ,(t)\ 



(2) 



A. Hamiltonian 

To study a CU2O3 ladder, we use the following three- 
band extended Peierls-Hubbard Hamiltonian, which in- 
cludes both electron-electron and electron-lattice inter- 
actions [l5llT^|: 



<ij>cr 



We impose periodic (open) boundary condition in the x 
(y) direction, i.e., there are two periodic Cu-0 chains 
along the x-direction (we term this oxygen O x subse- 
quently), connected together through the other oxygens 
(O y ). In this Hamiltonian, c\ a creates a hole with spin 
a on site i, and each site has one orbital (d x 2_ y 2 on Cu, 
and pa; or p y on O). The Cu (O) site electronic energy 
is €d (e p ). Uij represents the on-site Cu (O) Coulomb, 
Ud (Up), or the nearest-neighbor repulsion, U p d- The 
electron-lattice interaction modifies the Cu-0 hopping 
strength linearly through the oxygen displacement Uij: 
t P d(uij) = t pd (l ± oiUij), where +(— ) applies if the Cu-0 
bond shrinks (stretches) for a positive Uij ; it also affects 
the Cu on-site energies Cd(uij) = ed + PJ2j where 
the sum is over the three neighboring O ions. Other oxy- 
gen modes (buckling, bending, etc) are assumed to couple 
to electron charge more weakly and are neglected here for 
simplicity, but can be included as necessary within the 
same approach. We use variations around the follo wing 
set of model parameters used in 2D Cu02 models [lall4j : 
e p - e d = 4 eV, U d = 8 eV, U p = 3 eV, U pd = 1 eV, and 
K = 32 eV/A 2 , a = 2.0 eV/A, = 1 eV/A; we vary 
t p d = 1 ~ 5eV to simulate the pressure effect. This 
is clearly an oversimplified representation of pressure ef- 
fects, but serves to demonstrate the ground state phases 
and transitions. Effects of varying the coupling strength 
are also considered below with similar results. To ap- 
proximately solve the above model, we use unrestricted 
Hartree-Fock combined with an inhomogeneous general- 
ized RPA to study linear fluctuations of lattice, spin or 
charge m a supercell of size N x X 2 (we take N x = 8 
here) . 



B. Phonon spectral function 

The output of the calculation is the Hartree-Fock 
ground state and the linearized fluctuation eigen- 
frequencies and eigen- vectors with respect to it. From 



where e labels the five ions in the unit cell of the ladder: 
(1) O x ions in the lower leg, (2) O y ions in rungs, (3) O x 
ions in the upper leg, (4) Cu ions in the lower leg, and (5) 
Cu ions in the upper leg; Here the position of each ion is 
expressed by rj (t) — xi + d e + u e t (t) , where each of the 
terms represents the location of the Z-th unit cell origin 
Xj(= xil x ), time-dependent vibrational component u|(t), 
and position within the unit cell d £ : 



dW=6, d( 5 )= a L. 



(3) 



As noted above, for simplicity we consider Cu ions as 
fixed, and the motion of O ions oriented along the cor- 
responding Cu-0 bond: = u\e^ with ei = §3 = l x , 
62 = l y , §4 = §5 = 0. The scalar displacements can 
now be expressed in terms of the normal modes z n as 
u\(t) = ^Z n ^% un z n (t). Making a second-order expan- 
sion in the oxygen displacements, we obtain 



5(k,w) = 

£{[*2I*£»I 



3 i*£j a + 



+ [k x k v e -><*--*»>f + cc 



+ 



2muj, 



-[(I + n B )8(uj - u n ) + n B S(uj + w„)].(4) 



Here, <I>f 



,/T 



1 



<K,n ~ . ^xj.n: and n B = (e" 

is the thermal population of the phonon mode n. This 
is a generalization of the usual neutron scattering inten- 
sity expression [17j for the case of phonons with a larger 
real space unit cell. We plot 5(k, cj)/|k| 2 for k-directions 
sampling longitudinal modes, consistent with the com- 
mon experimental convention. 



C. Electron spectral function 

To investigate the neutral electronic excitations, we 
calculate the spectral function [il| : 



I (*o|0(k)|¥„) | 2 5{oj - E n - E ), 



(5) 



where |\& ) (\^f n )) is the RPA ground (excited) state 
whose energy is represented by Eg (E n ), and O(k) is 
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an operator, e.g. spin S(k) or charge n(k), summed over 
Cu- and O-sites: 

5 

O(k) = ^O e (k) e - lkde . (6) 

e = l 

The effect of an infinitesimal external field corresponding 
to the excitation ^f n can be represented by the change of 
an observable (O) in the state * = * + V^n (M < 1) : 

(O)~(O) + 5(O) n , (7) 

S(0) n CX (* |O|*n), (8) 

where (O)o is the expectation value with respect to the 
ground state. 

As noted earlier, we simulate the effect of pressure by 
varying t p d- We identify a transition from the SP to 
the BP state with distinct electronic and phononic sig- 
natures. Note that the BP states here are polaron bound 
states on different legs of the ladder (in contrast to same- 
chain BPs). 

III. RESULTS 
A. Polaronic ground states and phase transitions 

We first show in Fig.[T]the configurations of the ground 
states obtained by the Hartree-Fock calculation for sev- 
eral values of t p d- For 1.0 < t p d < 5.0, there are six types 
of groundstate configurations. (A) in Fig. ^ is a single- 
polaron (SP) state, which has n polarons in a staggered 
arrangement for an n-hole doped system (n = 4 here). 
(B) and (C) are diagonal bi-polaron (DBP) states, which 
have diagonally-bound polarons; and (D) and (E) are 
vertical bi-polaron (VBP) states, which have the same 
structure as a short segment of a vertical stripe in the 
2D system. (D) consists of site-centered VBPs, and (E) 
of bond-centered VBPs. (F) is the uniform (UNI) state, 
which does not have any local spin or charge modulation, 
or lattice displacement. 

Note that the undoped system shows the AF configura- 
tion without any lattice displacement (not shown here). 
One can expect that a singlet solid is the more likely 
exact undoped ground state, however the Hartree-Fock 
calculation favors an AF. 

With £7(a), ^(b)j — , ^(F) the energies of (A), (B), 
(F), respectively, we compare these energies, and de- 
termine the regions of t p d corresponding to the (A)-(F) 
phases for 1.0 < t p d < 5.0. The ^-dependence of these 
groundstate configurations is shown in Fig. [21 

There are five groundstate transitions for 1.0 < t p d < 
5.0. These are first-order transitions except for the one 
between (E) and (F), which is of second order. The main 
feature here is that the larger t p d, the more delocalized 
the ground state becomes. This follows from the fact that 
the transition with increasing t p d is SP — > BP — > UNI. 
By studying the polaron eigen-functions, we find that 
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FIG. 1: (Color online) ground states for several values of t p d'- 
(A) t pd = 1.2, (B) t pd = 1.6, (C) tpd = 1.8, (D) t pd = 3.0, 
(E) tpd = 3.5, and (F) t pd = 4.7. The circles (radius) and the 
arrows (length and direction) represent the hole and spin den- 
sities, respectively. Lines at the positions of O-ions represent 
static displacements, and shadows show where the polarons 
are located. 



the transitions occur when the SP-SP or BP-BP overlap 
achieves sufficient levels, resulting in the polaron melting 
,i.e., the sequence of the transitions SP — > BP — > UNI. 

The transition DBP — > VBP with increasing t p d is 
similar to the transition from diagonal to vertical stripe 
with increasing hole-concentration in some of the high-T c 
cuprates with Cu02 planes 0,|2(J. The energy differ- 
ence between (D) and (E) near the transition point is 
smaller than for the other transitions. -E'(d) — -^(E) an d 
E(e)—E(f) both asymptotically vanish, as t p d approaches 
the (E)-(F) transition point from smaller t p d- This tran- 
sition is associated with recovery of the broken symmetry. 
Below, we show the existence of the VBP sliding mode, 
which recovers the translational symmetry. 

B. Sliding mode in bi-polaronic states 

The ipd-dependence of the phonon eigen-frequencies 
obtained by a direct-space RPA calculation is shown in 
Fig. |21 A main branch lies in the range of 80-85 meV 
and is consistent with the main phonon branch observed 
experimentally. This frequency range is insensitive to the 
doping level. However, some characteristic local modes 
are induced below this frequency range by hole doping, 
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FIG. 2: fpd-dependence of the energy difference between the 
(meta-) stable states with a = 2.0: -E(a) — -E(es) (filled circles), 



J (B) ■ 



■E(C) (open squares), £7c) — ^(d) (open circles), 



-E(e) (filled squares) , and S( E ) — S( F ) (open triangles). E, 



Ho)' 

(D)- 



i?(E) is of order 10 



3 meV per unit cell (CU2O3), while the 
others are of order 10 _1 ~ 1 meV. Each state (A)-(F) is the 
ground state in the corresponding region separated by the 
thick grey lines. The lines are guides to the eye only. 
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FIG. 4: (Color online) (a) charge excitation for k= (-§,0) at 
t p d = 3.5 eV from the RPA calculation, (b) Excited state at 
E ~ 13 meV. The ground state is (E) in Fig. Q Comparing 
this excited state with the ground state, a sliding of BPs is 
found, (c) Energies of soft phonon modes (open circles) and 
charge sliding modes (filled circles). 



similar to the 2D cases [H E3. Additionally, an ex- 
tremely soft mode is found in the (D) and (E) phases. 
The frequency of the soft phonon mode in the (E) config- 
uration is purely imaginary in the (D) phase region, and 
vice versa. To further understand this soft mode, we cal- 
culate the corresponding electronic excitations [Fig. Q]. 
We find one soft charge excitation [Fig. 0] (a)], which 
shows a sliding mode [Fig.@|(b)] and whose excitation en- 
ergy shows the same behavior as that of the soft phonon 
mode. We identify the soft phonon mode as one cou- 
pling with a sliding mode of VBPs along the ladder. We 
identify the frequency of this mode as a pinning energy 
of the VBP sliding, corresponding to a Peierls-Nabarro 
barrier from the lattice discreteness. In the (D) phase 
[(E) phase], VBPs are pinned Cu site-centered [bond- 
centered] by a potential energy; the potential energy is 
minimal at a Cu site [an 0^ site] and maximal at an O x 
site [a Cu site]. The minimum and maximum points of 
the potential energy exchange at the (D)-(E) transition 
point. We will discuss this further below. 

From the softening of the sliding mode and the second- 
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FIG. 3: ^-dependence of the phonon eigen-frequencies. The 
lines are guides to the eye only. 



order-like behavior of the energy differences, -E^d) — ^(E) 
and E(e) — En?) near the (E)-(F) transition point, we 
conclude that the pinning potential becomes flat at that 
transition point. 



C. Effects of electron- lattice coupling 

We next consider effects of the electron-lattice cou- 
pling on the groundstate configuration. For this purpose, 
changing the coupling strength a (=1.0, 2.0, 3.0, 4.5) we 
calculate the critical values of t p d in the same manner as 
above. In this way, the t pc i — a phase diagram is found as 
in Fig.Q] Increasing the coupling strength tends to raise 
the critical values of t pc i except for that between (E) and 
(F): The boundary between (E) and (F) is insensitive to 
the change of the coupling strength. The difference in the 
a-dependence shows that the transition between (E) and 
(F) has a different character than other transitions: As 
seen above, the transition between (E) and (F) is most 
likely of second-order. 
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FIG. 5: phase diagram obtained from the calculation of 
the critical values of t v d for various electron-lattice coupling 
strength a. (A)-(F) correspond to the states in Fig. Q The 
lines are guides to the eye only. 
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IV. DISCUSSION 



We now discuss the details of the results obtained 
above. First, we analyze the polaron size by using Gaus- 
sian fitting, and describe the overall picture of the tran- 
sitions as the derealization of polarons induced by the 
pressure. Next, we discuss the transition between site- 
and bond-centered VBPs and the pinning frequencies of 
a CDW in terms of Ginzburg-Landau theory. 



A. Size of polaron 




(E)a=4.5,t pd =4-0 
1 8 "1 2 3 4 5 6 

Cu site i 

FIG. 6: hole-density profile fitted by Gaussian. The left pan- 
els for a = 2.0, and the right for a — 4.5. The upper, middle, 
and lower panels for (A), (C), and (E) states, respectively. 
Circles show the hole-density averaged over nearest O sites 
using the data of the the Hartree-Fock calculation. Here the 
profile along only one of the ladder legs is shown. 



Fig.|Hlshows hole-density profiles for several cases. The 
data shown by circles in Fig. includes averages over 
nearest O sites: 



E 

j: nearest 
neighbors 



noj, 



(9) 



where nc u ,i and n-oj are hole densities, respectively, at 
Cu and O sites obtained by Hartree-Fock calculation. We 
fit the data by a least-squares method with the following 
function: 



(x 



2a 2 



(10) 



Here, the center of the i-th polaron is represented by 
{xi,yi), and h takes a value of about 1.1. In the SP 
and VBP cases (upper and lower in Fig. EJ), the data is 
well-fitted by the function in Eq. (|1U|I . In the DBP case, 
on the other hand, it is not as well-fitted, especially for 
large a (middle right in Fig. . We define the size of 
a polaron as 2a /a. The t p d dependence of this polaron 
size is shown in Fig. {7\ In the (A) case, we find 2a /a 



D 2.5 



(D),(E) 




FIG. 7: tpd dependence of polaron size at a — 2.0. (A) 
filled diamonds, (B) open triangles, (C) open circles, (D) open 
squares (E) filled triangles. (B) and (C), and (D) and (E) lie 
almost on the same line. 



is between 1.1 and 1.2 within the t p d range between 1.0 
and 2.5. There is little change with increasing t p d in this 
case. The (B) state shows the same behavior as the (C) 
state. The polaron size of the DBP states increases as t p d 
increases. This polaron size growth in the DBP states is 
greater than in the SP state. However, the polaron sizes 
in the VBP are larger than those of SP or DBP. 

Comparing with Fig. |3 we find that the phase transi- 
tions occur from a state of smaller polarons to another 
state of larger ones. Especially, we can understand the 
aspects of the phase transitions (D)— >(E)— >(F) by con- 
sidering the pinning potential. Derealization of the po- 
laron induced by changing t p d causes less distortion of the 
lattice, as shown in Fig. 2] Therefore, the potential en- 
ergies for site- and bond-centered states should change: 
For small t pc t, the pinning at Cu is stronger than that 
at O. As is discussed in Appendix ^] if the change of 
t p d varies the ratio of the pinning potential at Cu and 
O sites, then phase transitions are induced. The phase 
transitions are also well-described by Ginzburg-Landau 
theory. This is discussed in the following section, where, 
within the Ginzburg-Landau picture, we explain the be- 
havior of the pinning frequencies. 



B. Pinning frequency 

We now discuss the transition between site- and bond- 
centered VBPs and the behavior of the pinning frequen- 
cies in terms of a Ginzburg-Landau theory. We can de- 
scribe the aspects of the transitions between VBP and 
UNI, and between site- and bond-centered VBPs by in- 
troducing the Landau function (See Appendix iBl for de- 
tails): 



F = dx f[ip(x), m{x)\. 



(11) 



Here i(>(x) and m(x) are the charge and spin order pa- 
rameters. m{x) is the staggered spin density, and ip(x) is 
defined by the deviation of the charge density from the 
uniform state: 



m{x) 
ip(x) 



(-l)iS(x) 

p(x) ~ po- 



rn) 

(13) 



In general, / can be written in the following form for an 
L-site-periodic commensurate CDW plf : 

f[ip(x),m(x)} = f [ip(x),m(x)] 

+ P {x)iP(x) L + q(x)^(x) 2L . (14) 

Here, /o is concerned with the lattice-independent spin 
and charge ordering, and the remaining terms lead to the 
lattice pinning effect of the CDW. We assume fo is of the 
form 

fo[ip(x), m{x)] = rom(x) 2 + uom(x) 4 

+ s m(x) ip(x) + v ifj(x) 2 (15) 

with vq > 0, that is, the charge order is induced by the 
magnetic order |22j | . If we write the charge order param- 
eter in sinusoidal form, the amplitude p\ is found to be 
(see Appendix |B|) 



Pi oc t — t p d, 



(16) 



where to is the BP-UNI transition point. 

Considering the small oscillations around the equilib- 
rium state, the pinning frequencies are derived in Ap- 
pendix [5] as 



\t 



pd 

|io-i P d| (i - 2)/2 fort pd <t 



y 1/2 



for t 



pd 



tr 



(17) 



where t c is the transition point between site- and bond- 
centered VBPs. Both expressions are plotted in Fig. 
(b), and show a good agreement with the RPA data. If 
t c is not very far from to, the following form well-describes 
the behavior of the phonon frequencies as a function of 
t p( i over the whole region around the transition points 
[see Fig. El (a)]: 



oc \t p 



t, 



a/2 



|*0 — tpd\ 



(L-2)/2 



(18) 



As shown in Appendix [Bj the t p d dependence of the 
energy difference between site- and bond-centered VBPs 
is given by 



\AF\ 



tpd 
1*0 



tpd 



for tpd 
L for tpd 



tc 

to 



(19) 



These functions are plotted in Fig. 0(c), where \Em\ — 
E(e) I from the Hartree-Fock calculation are also plotted 
for comparison. 

Eq. (|17|l well reproduces the features of the pinned 
CDW in the L = 4 case. Here we investigate energies and 
the pinning frequencies for different dopings (resulting in 
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FIG. 8: Comparison of the functions obtained from Ginzburg- 
Landau theory with the results from Hartree-Fock and RPA 
calculations for the i-hole doping (L = 4) case. Lines show 
the functions (the explicit forms are shown below). The fre- 
quencies from the RPA calculation are plotted with open tri- 
angles for tpd^tc, open circles for t p d?L,t c , and filled circles for 
tpdr^to- The energy difference \Er D \ —E(e)\ from the Hartree- 
Fock calculation are plotted with open diamonds for t p d^t c , 
open squares for t v d^jt c , and filled squares for t p d^to- (a) Pin- 



ning frequencies: S7 oc \t p d — t c \ \to — t p d\ ~ ■ (b) Pinning 
frequencies (log-log plot): SI oc (to — t p d)^ L ~ 2 ^ 2 for £ p( j~io 
(solid line) and Q, tx (t c — tpd) 1 ^ 2 for t p d ~ t c (dashed line), 
(c) Energy difference between site- and bond-centered VBPs: 
|AF| oc (to—t p d) L iortpd^to (solid lines) and |AF| oc (t c —t p d) 
for tpd ~ t c (dashed lines). 



different-period CDWs), and further show the validity of 
Eq. for other L. Figs. M and show the results 
for the i-hole doping (L — 3), and the ^-hole doping 
(L = 5) cases. In the both cases, the results suggest that 
Eq. (|17fl agrees with the RPA calculation: The Hartree- 
Fock and RPA calculations were performed in the sys- 
tems whose sizes were 6 x 2 (L = 3) and 10 x 2 (L = 5). 
The tpd dependence of the energy difference between site- 
and bond-centered VBPs is also well-described by the 
function obtained from Ginzburg-Landau theory in both 
cases. 

Similar to the i-hole doping, for i-hole doping. 
Eq. i|18|) is also a good approximation. For |-hole dop- 
ing, Eq. Ijl 80 does not give a good agreement with RPA 
data. This is not surprising, since t c is far from to in this 
case. 



V. SUMMARY 

In summary, we have modeled a pressure effect in 
a CU2O3 ladder system by using a multi-band Pcicrls- 
Hubbard model and simulating the effect of pressure 
through the hopping strength t p d- With increasing t p d, 
we find a sequence of transitions from SP/BP charge lo- 
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(a) (b) 




FIG. 9: Pinning frequencies and energy difference for the |- 
hole doping (L = 3) case. See the figure caption in Fig. |H] 



(a) (b) 




FIG. 10: Pinning frequencies and energy difference for the 
i-hole doping (L = 5) case, (a) fi = C(to - t vd ) [L ' 2)/2 for 
t p d around to (solid line), D. = C\t c — tpdl 1 ^ 2 for t p d around 
t c (dashed line), and the RPA data of the pinning frequencies 
(same as in Fig. [SJ . (b) and (c) are the same as in Fig. |S| 



between site- and bond-centered VBP states, although 
the energies are very close. Calculations of the phonon 
eigen-frequency and electronic excitation in the VBP 
phase yields a sliding mode of VBPs with weak pinning. 
The pinning energy in the bond-centered VBP phase is 
around 15meV at most for t pc t ~ 3.5 eV. Increasing t p d 
up to t p d ~ 4.6 eV makes the pinning zero, and a tran- 
sition from the VBP to UNI state occurs. These results 
suggest experimentally exploring pressure de pen dence of 
the low-energy modes found by Kitano et al. [T3, EH and 
Blumberg et al. 01 ^ these modes correspond to those 
we have identified, then their pressure-dependence will 
follow the interesting pattern in Fig. |3] Studying IR, 
Raman and optical signatures would further clarify the 
mode assignments. The resonant soft X-ray scattering 
technique of Refs. |8|, |9| could also be used to probe our 
predicted charge ordering structures as a function of pres- 
sure. 

The sequence of phases (SP-BP-UNI) is reminiscent 
of the insulator-metal transition with doping observed 
in other low-dimensional broken-symmetry groundstate 
materials, including conjugated polymers [23j and lay- 
ered cuprates. It is tempting to associate the mode soft- 
ening with the onset of a sliding CDW in the spirit of 
Frohlich. However, as the phonon-fluctuations soften, 
additional degrees of freedom (quantum lattice and spin 
fluctuations) become relevant and need to be considered 
— in particular to identify the superconductivity mecha- 
nism. The superconductivity is observed experimentally 
in a finite range of pressure. Whether this can be asso- 
ciated with the finite range of t p d with low pinning fre- 
quencies (Fig. |3J) requires comparison with more detailed 
experiments, but is a tempting scenario. We have also 
explored other values of a and found the same general 
phase sequence as a function of t p d shown in Fig. [SJ 

In this study, we mainly considered only the ^-hole- 
doped case and related commensurate dopings. Other 
cases including incommensurate fillings with discommen- 
surations will be reported elsewhere. We can expect 
related transitions with doping as with t p d, since they 
should both be controlled by SP or BP wave-function 
overlaps. 



calization to sliding to delocalization, all occurring within 
a magnetically ordered background. The ground state 
has the same number of SPs as doped holes in the case 
of small t p d (t p d < 1.6 with the parameters used here). 
The ground state configuration changes as SP — * BP — > 
UNI states, as t p d increases. While SPs are localized and 
isolated, BPs are partially delocalized. This means the 
pressure produces a more delocalized ground state. In the 
BP phase, there is also a phase transition between DBP 
and VBP states. A similar transition has been found 
in some other cuprates: namely, the transition between 
diagonal- and vertical-stripe states induced by hole dop- 
ing. 

In the VBP phase, we also find a soft mode transition 
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APPENDIX A: PINNING POTENTIAL 

Here we attempt to describe the pinning potential of 
VBP states, and show the details of the discussion re- 
garding the transition between site- and bond-centered 
VBPs in Sec. IIV Al Since we are interested only in VBP 
states with different phases here, the system we consider 
can be reduced to one dimension. Therefore, the pin- 
ning potential is a function of x, and has a minimum at 
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x = na for the site-centered VBP and at x — (n + |)a 
for the bond-centered VBP. Approximately, the pinning 
potential may be attributed to the potential at Cu and 
O sites: 

E(x) = (£cu,i(z) + £o,i(x)) , (Al) 

i 

where the phase is chosen such that the state becomes 
site-centered when x — 0. We continue our discussion 
with the following two assumptions: (1) both £qu and 
So are composed of Gaussians, 




and (2) both a\ and 02 are comparable to a. The case 
we consider here can satisfy these conditions. Since the 
density profile of VBP states is well-fitted by Gaussians 
(shown in Sec. IIV|I . we expect the contribution of the 
partial free energy of the Cu and O sites to the total is 
also formed of Gaussians, and that the size of the Gaus- 
sians should also be similar to the VBP size. 

We investigate the potential energy in Eq. (|Aip with 
assumption (1). The pinning potential can be expanded 
as a Fourier series: 

E(x) = J2 E n cos(nGx) , ( A4) 

n 

where G = — . The Fourier coefficients are given by the 
a 

following form: 



E n — — 



2"KG\ 



27T(7l 

dx exp 



exp 



2af~ 



+ 



C 2 



27TCT; 



■ exp 



(x - ai+ f ) 2 



2a 2 



cos (nGx) 



x 



., ( cos(nGx) + (-1)' 



Co 



27T(72 



dx exp 



2a\ 



cos (nGx) 



"Cx I 


' n 2 G 2 a 2 \ 


CXp<! 









+ (-!)" — exp {- 



n 2 G 2 a 2 



(A5) 

(A6) 
(A7) 



From Eq. I)A7|1 . it follows that E n for large n vanishes. By 
considering the fact that the Gaussian almost vanishes 
at three-fold half maximum full-width, we evaluate the 
condition to neglect the components as: 



n 2 G 2 a 2 2n 2 n 2 a 2 3 2 



, 2 > 2 



(A8) 



Here all variables are positive, so that the condition is 
given by 



2a 1 
— > -. 

a 11 



(A9) 



If we consider the case that both o~\ and o~2 are compa- 
rable to a (assumption (2)), E(x) is approximately rep- 
resented by the cosine curve or a slightly modified one: 

E{x) - E + E 1 cos (Gx) + E 2 cos (2Gx). (A10) 

Here E2 is small and it does not change the shape of 
cos (Gx) very much unless E\ is small as well as E2 ■ In 
such a situation, only one of the site- or bond-centered 
VBP states is stable, and is determined by the relation 
between the C\ and C2 magnitudes. However, if E\ is 
smaller than AE2, E(x) has minima at both Cu and O 
sites. This can happen when C\ ~ C 2 (see Fig. ITT1 (a,)). 

In Fig.^](b), the scheme of phase transition which fol- 
lows from the Hartree-Fock calculation (Fig. [2J) is shown. 



From the phonon frequency calculation (Fig.[3J), the bi- 
stability interval, t\ — t%, is very narrow and difficult 
to identify numerically. Comparing this scheme and the 
behavior of dE(x) / dx, we understand the phase transi- 
tion (D)— >(E)— >(F) as follows. (1) Since the site-centered 
state (D) is found at small t p d (< ^2), C2/C1 has to be 
less than 1 in this case. (2) Since the double-minimum 
region is very narrow (t\ ~ £2), the zig-zag boundary near 
C2/C1 ~ 1 in Fig. ^2 is almost flat when t p d is close to 
the transition point. (3) Since the bond-centered state 
(E) is found above the transition point (> ti), C2/C1 
should be found in the upper region (> 1) across the co- 
existence point t c . (4) Increasing t pc i far above ti, E(x) 
becomes close to a constant Fq, which implies that the 
transition from (E) to (F) occurs. 



APPENDIX B: GINZBURG-LANDAU THEORY 
OF THE SITE-CENTERED TO 
BOND-CENTERED VBP TRANSITION 

Here we give details of the discussed transition between 
site- and bond-centered VBPs in Sec. IIVBI First, we 
describe the statics of the transition by considering the 
free energy. Then, by considering small long-wavelength 
oscillations around the ground state, we estimate the pin- 
ning frequencies. 
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(a) 




(b) 



'c 



O Cu O Cu O Cu O 



FIG. 11: (a) The derivative of E(x). — and + show sign of 
dE(x)/dx. Blackened area represents dE(x)/dx ~ 0. Here 
both (Ji and 02 are taken as 0.4. If we take larger a\ and 02, 
the zig-zag boundary becomes flatter, (b) Scheme of transi- 
tions in bi-polaronic states. Solid lines show the center of bi- 
polarons at each t p d. In the case t p d < ti (ii < t p d < to), the 
system shows the site-centered (bond-centered) VBP ground 
state. For £2 < t p d < ti, the pinning potential has two min- 
imal points, and at t c there is coexistence of site- and bond- 
centered VBP states. In to < t p< j (blackened zone), the (F) 
state becomes the ground state. The larger ai and 02, the 
narrower the double-minimum zone (ti < t p d < £2) becomes. 



We start with the Landau function: 



F = I dx [ rom(x) 2 + uom(x) 4 



+ Som(x) ip(x) + v ip(x) 

+ p(x)ip(x) L +q(x)<ip(x) 2L ] (Bl) 

We write p(x) and q(x) in an expanded form, for example 

00 

p(x) = J^Pn cos(nGx). (B2) 



As seen in Sec. IIVI the density profile of the VBP is well 
fitted by large-size Gaussians, so we can neglect higher 
harmonic terms of the order parameters and describe 
them by sinusoidal waves: 



ifj{x,(j)) ~ pi cos(^x + |), 
m(x, 4>) ~ mi sm(j^x + ■$£). 



(B3) 
(B4) 



We may constrain p\ > 0. 

We evaluate the free energy using the following integral 
(I + 0): 



f ^- co ^(f x +l) C os(lGx) 
Jv V 



(n-l-5 n )/2 



E 

h=0 



1 [n 



2 n \h 



S n -2h,lL cos(l(/)), (B5) 



where 



6 is Kronecker's delta, and 5 n is 



a/ {n-l)\l\ 

unity for even n and zero for odd n. Integrating over 



volume V, the free energy per unit volume is 
Ffiitpd, p\,m\) = F (0) (i p(i ,pi,mi) 

+ [PiPx + qip\ L ] cos(0) 
+ g 2 p? i cos(20), (B6) 
F(°\t pd ,pi,mi) = f Q ml+u mj 



s m\pi + vqpI, (B7) 



where 



Pi 



u 



2LP1, 9i 



9i = 25rgi, 



gM , v = 2 u o, r = 2^0, «o = 4 



92=2^92, (B8) 

ho, (B9) 



Note that the signs of these variables are the same with 
or without tilde. The ro and uq terms govern the BP- 
UNI continuous transition, and the rest characterize the 
transition between site- and bond-centered VBPs. Sub- 
sequently, we neglect the qip 2L term in Ea. (|B6|) . since 
this term would be smaller than the p\p\ term. 

First, we consider the BP-UNI transition. Phenomeno- 
logically assuming r oc t p d — io (to is the BP-UNI 
transition point, where p\ — 0) and tto > 0, we find 
Tfi\ oc (to — tpd) 1 / 2 for t p d < to. From the last two terms 
of Eq. CEZI), it follows that 



pi cx m 1 oc t - tpd 



(B10) 



for tpd < to- Next, we consider the transition between 
site- and bond-centered VBPs. In the case t p d ^ t c (t c 
is the coexistence point), the minimum point of this free 
energy is controlled by the p\ term. At t p d — t c , the 
Pi term vanishes, and the minimum of the free energy is 
determined by the qi term. 

The minimum point is given by 



7r for p\ > 
for pi < 



(Bll) 



<f> = tt corresponds to the bond-centered case at t p d > t c , 
and <j) — to the site-centered case at t p d < t c . At the 
critical point (t p d = t c ), the p\ term vanishes; and for 



q2p\ L < 



(B12) 



there are minima at both = and 7r. 

Next we consider the small oscillations around the 
equilibrium state: cf> — 4>o + S<p (4>o takes either or 
7r for the ground state). Using pi cos(0o) = — \pi\ an d 
cos(2</)o) = 1, the free energy is expanded for S(f> as 



F4, 0+ 84,(tpd,tp,m) 

= F°(tp d ,pi,mi) +P1P1 cos(^o) 

(254>r 



q 2 p\ L cos(20o) 

[F°(t pd ,pi,mi) - 
1 ,, 



1 - 



\Pi\Pl + 



hp\ L \ 



(B13) 



(B14) 
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The second line of Eq. (|B14|) is used to find mode fre- 
quencies. The Lagrangian for the oscillation of VBPs is 
given by 



C = 



M d{84>) MQ 2 



m 7 



const. 



2 dt 2 
Here, M is the effective mass of the CDW, and 



(B15) 



(B16) 



Using pi oc t p d — t c and p\ oc to — t p d, and supposing 
M oc p\ , for tpd < to, we find 



/ i j tpd tc | ^ for tpd ^ t c 

\to-t pd \ (L - 2)/2 ioxt pd <to 



(B17) 



It also follows from Eq. (|B6|I that the t p d dependence 
of the energy difference between site- and bond-centered 
VBPs, |AF|, is given by 



"■- ''^-'^^-t^t'C-t- (B18 » 
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